arXiv:1509.02126vl [hep-th] 7 Sep 2015 


CERN-PH-TH/2015-182 


The spectrum of anomalous magnetohydrodynamics 


Massimo Giovannini 


0 


Department of Physics, Theory Division, CERN, 1211 Geneva 23, Switzerland 
INFN, Section of Milan-Bicocca, 20126 Milan, Italy 


Abstract 

The equations of anomalous magnetohydrodynamics describe an Abelian plasma where 
conduction and chiral currents are simultaneously present and constrained by the second law 
of thermodynamics. At high frequencies the magnetic currents play the leading role and the 
spectrum is dominated by two-fluid effects. The system behaves instead as a single fluid in 
the low-frequency regime where the vortical currents induce potentially large hypermagnetic 
helds. After deriving the physical solutions of the generalized Appleton-Hartree equation, 
the corresponding dispersion relations are scrutinized and compared with the results valid 
for cold plasmas. Hypermagnetic knots and fluid vortices can be concurrently present at 
very low frequencies and suggest a qualitatively different dynamics of the hydromagnetic 
nonlinearities. 


^Electronic address: massimo.giovannini@cern.ch 



1 Introduction 


Electrically conducting media are customarily described as a single fluid in the low-frequency 
branch of the plasma spectrum. This approach has been extensively applied to the analysis of 
hydromagnetic nonlinearities [T] evolving in terrestrial |2] and astrophysical plasmas mw- 
The same strategy cannot be extended to higher frequencies where the one-fluid description 
is no longer tenable [6] and the plasma must be treated, at least, as a double fluid. This well 
known aspect of conventional electromagnetic plasmas stems directly from the properties of 
the vector currents which are associated, in the high-frequency limit, with the ions and with 
the electrons. When the plasma is globally neutral the total vector current is instead Ohmic 
in the low-frequency domain. 

A problem of similar nature occurs in anomalous magnetohydrodynamics describing 
a charged fluid where axial and vector currents are simultaneously present: while the axial 
currents are not conserved because of the triangle anomaly, the vector currents are even¬ 
tually Ohmic. The purpose of this investigation is a systematic discussion of the spectrum 
of anomalous magnetohydrodynamics (AMHD in what follows). The equations of AMHD 
differ from the ones where only chiral currents are present IHli at finite fermionic density. 
They generalize the system firstly explored in Ref. [10] accounting for the evolution of the 
hypermagnetic and hyperelectric fields in the electroweak plasma. Indeed, in the symmetric 
phase of the electroweak theory the non-screened vector modes of the plasma correspond 
to the hypercharge which has a chiral coupling to fermions. The axial currents may be 
associated with the evolution of the chemical potential or with the presence of an axionlike 
held [ill 112] (see also m)- In both cases the plasma may host parity-odd configurations 
of the gauge fields characterized non-vanishing hypermagnetic gyrotropy B ■ V x B which 
is the hydromagnetic analog of the kinetic gyrotropy (i.e. v -V x v) naturally appearing in 
the discussion of mean-held dynamos oa. The dynamical production of hypermagnetic 
knots and Chern-Simons waves during inflation offers a potentially viable mechanism for the 
generation of the baryon asymmetry of the Universe (see last two papers in [ID])- In AMHD 
the hypermagnetic current is complemented by a vortical current possibly leading to the 
formation of fluid vortices. 

The same class of physical systems previously discussed in the electroweak plasma also 
arises in the framework of the so called chiral magnetic effect [H] . Both phenomena are often 
presented as macroscopic manifestations of triangle anomalies. The model of chiral liquid 
emerging in the context of AMHD could then be relevant also in the context of the chiral 
magnetic effect insofar as axial currents and quark vector currents are concurrently present in 
the strongly interacting plasma. In the absence of finite conductivity effects (see e.g. mm) 
the validity of the second law of thermodynamics is guaranteed by the simultaneous presence 
of an hypermagnetic current and of a chiral vortical term. In AMHD the vector currents 
(eventually responsible of Ohmic dissipation), the chiral currents (determining the anomalous 
effects) and the vortical currents (required by the second principle of thermodynamics) are 
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all described by the appropriate kinetic coefficients. Whenever possible AMHD will be 
discussed in analogy with the spectrum of conventional plasmas. More specifically the plan 
of this investigation is the following. In section [2] we discuss the relativistic problem and 
derive the general form of the kinetic coefficients. In section |3] the properties of the two-fluid 
equations are analyzed while section 0] is devoted to the dispersion relations in the high- 
frequency domain. The one-fluid equations and their implications are presented in section 
[5l Section |6] contains our concluding remarks. To avoid digressions some relevant technical 
aspects have been relegated to the appendix. 

2 The relativistic formulation and the total entropy 

The conservation of the total energy momentum tensor and the evolution of the chiral and 
vector currents determine the relativistic form of the second law of thermodynamics. If the 
four-divergence of the entropy four-vector is to be positive semi-definite (as implied by the 
generalized second law) the chiral and vector currents must contain supplementary kinetic 
coefficients corresponding to the hypermagnetic and to the vortical currents. In what follows 
the dissipative effects are included in the framework of the Landau approach: the total four- 
velocity coincides then with the velocity of the energy transport defined from the mixed 
components of the total energy-momentum tensor. 

2.1 Ohmic and chiral currents 

In the simplest situation the total energy-momentum tensor of the system what 

follows) consists of four qualitatively different contributions: the energy-momentum ten¬ 
sor of the charged species (denoted by the energy-momentum tensor of the chiral 

species (labeled by T)(^)), the dissipative contribution and the gauge contribution 

(corresponding to an hypercharge gauge-field strength Y^y)\ 

rp{tot) _ rp{ + ) I rp(-) . rp(R) , rp(diss} , q- (9 

flU ~ 111/ ' IIU ' III/ ‘ Iiu ' 

The covariant conservation of implies, as usual, that = 0 where denotes 

the covariant derivativ^ defined from the metric tensor (with signature mostly minus). 
The chiral and the conduction currents coexist but are not bound to coincide: they obey 
different equations. More specihcally the anomalous current is not covariantly conserved and 
its evolution can be written a^ 

= = + ( 2 . 2 ) 

^The discussion will be conducted in a general relativistic formulation even if the spectrum of AMHD 
will be discussed in flat space-time. 

^Note that nn and are respectively the concentration and the four-velocity of the chiral species. 
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where Ar is a numerical factor that is determined by the specihc nature of the chiral species 
and by the coupling to the hypercharge field; note that in the Landau frame i'rU^ = 0. 
Conversely the conduction current is covariantly conserved, it is a source of the evolution 
equations of the gauge fields and it may even contain a dissipative contribution: 

V^F^'^ = 47r/, V^j^ = 0, (2.3) 

where j'^ = {q±n± u± + z/±). The dual held strength obeys, as usual, = 0. 

This is the approach followed in [7] which differs from other more conventional approaches 
[la M\ where the anomalous current and the conduction current are identihed. In the 
present approach the anomalous current is not directly the source of the evolution of the 
hypercharge. 

In a general relativistic description and its dual account for the evolution of the 
hyper charge £eld0 however the gauge held strength can be decomposed into the hyperelectric 
and hypermagnetic parts denoted, respectively, by and 


La/3 = SaUf} - SfiUa + EafipaU^B^", (2.4) 

where Ea^pa = \/—g ^a/3pa and Sa/Spa is the four-dimensional Levi-Civita symbol while g is 
the determinant of the metric tensor. The total four-velocity of the system follows from 

{p + p)u>^u'' = Y.\PA + P{a)] «fa) (2-5) 

a 

where w = {p + p) denotes the total enthalpy; the sum in Eq. (12.hh runs over all the species 
of the plasma, both charged and chiral. Close to an equilibrium situation the four-velocity of 
the anomalous species coincides with the bulk velocity of the plasma and, therefore, ~ u^. 
The vorticity four vector can then be dehned as: 

( 2 . 6 ) 

From Eqs. (12.dh and (12. 6 h it follows that the four-divergences of and uj^ are given by: 

w'WpOJ^ = —2uj‘^daP — 2n£°'Ua, (2.7) 

wWpBf^ = 2wYp„uPu‘^ + Updo,pY^^ YUpY^pf?'^^, ( 2 . 8 ) 

wVpS^ = w[A7irua-Yf^PupUp] + Yf^"<upd^p + Y^'^upY^o.r, (2.9) 

where w = p + p is the enthalpy density of the huid. Equations fl2.7p . fl2.8p and fl2.9p have 
been obtained in the globally neutral case where h = h+ = h_ and q+ = q = —g_ but they 
can be easily generalized to the case where the plasma is not globally neutral. 

^As soon as we speak of hyperelectric and hypermagnetic fields we are implicitly assuming that the 
plasma has a finite conduction current so that a preferred frame can be selected where the electric fields are 
suppressed. Even if the electric and magnetic fields are non-relativistic concepts, it is practical to introduce 
the electric and the magnetic components of the gauge field strength in a generally covariant language. 
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2.2 First and second principles of thermodynamics 

Denoting with the chemical potential associated with the anomalous species, the hrst 
principle of thermodynamics demands: 

dE = TdS — pdV + fiiidNR, w = p + p = Tq + (2-10) 

The fundamental identity E = TS — pV + PrNr can be divided by a hducial volume 
and the result is the one reported in the second relation of Eq. fl2.10l) where is the 
entropy density and p the total energy density of the system. Combining the two relations 
of Eq. fl2.1Up further thermodynamic relations can be obtainecj^. Since the anomaly-induced 
currents are protected by topology they are not associated with dissipative effects. Thus the 
entropy production of the plasma must only come, in the relativistic case, from the viscosity 
coefficients or from the Ohmic contributions but neither from the chiral currents nor from 
the corresponding diffusive contribution. The absence of dissipative contributions stemming 
from the anomalous sector demands that the total entropy four-vector must be supplemented 
by two further coefficients S^, and Sb'- 

= + + SbB^. ( 2 . 11 ) 

The covariant conservation of the total energy momentum tensor can be written as 

rjnfliy 

= Z, (2.12) 

where we assumed, for the sake of simplicity, a global charge neutrality of the plasma and 
a corresponding Ohmic form for the charged species, namely = aY°‘'^ Ui, where P" = 

— u^u°‘ is the standard projector. The function Z appearing in Eq. fl2.12p is given by 

Z = V^^S^uj^ + SBB^^-'^Y^p-dpJIj,u^-ABjlByapY-^. (2.13) 

We remark that the specihc dehnition of the entropy four-vector depends on the chemical 
potential of the system. However, since the coefficient Ar does not have a dehnite sign, the 
anomalous currents may even lead to violation of the second principle of thermodynamics 
unless Z vanishes identically. 

2.3 Magnetic and vortical coefficients 

The vortical and the magnetic currents modify also the diffusive contributions denoted, 
respectively, by and in Eq. (12.131) . Four different coefficients parametrize the relation 
between (z/", z/)j) and (ca", H"): 

z/“ = + = Ar^u^ + Arb B'^, (2.14) 

®Like, for instance, <;daT + nudaRR = daP or daP = Tda'i + RndaUR- 
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where (A^^, A^) and Ajib) ah depend on the chemical potential and on the temper¬ 

ature. Using Eqs. fl2.7l) . fl2.8p and fl2.9p the condition Z = 0 together with the explicit 
expression of Z (see Eq. fl2.13p l becomes: 


25b- 

2 

- CTcUJ 

W 










- 1 - 


: 0 


(2.15) 


where Va and Qa are two differential operators defined respectively, as: 


Va = - —S^daP - daPn Aruj, Qa = daSs - —S^daP - daPR Arb- (2.16) 

w w 

The results of Eqs. (I2.13p - (l2.15p follow easily if we recall that, by dehnition, u^^Ua, u^S /3 
and u'^B^ are all vanishing. 

To satisfy the condition expressed by Eq. fl2.15p the four-vectors multiplying a;“ and 5" 
must vanish together with the coefficients of the terms multiplied by u°‘Ba and S°‘Bq. We 
then arrive at the following conditions: 


Va — 0, Qa — 0, Ar — 4:PrAr, A^ — 2TSb, 5aj — 0. (2-17) 


If, as established, 5^^ = 0 then Eq. fl2.15p also implies that A^uj = 0. All the coefficients we 
ought to determine depend on and on the pressure. Thus the conditions of Eq. fl2.17p 
are equivalent to the following system of equations: 

where A^^ = 2TSb and As = ‘^ArJI^. Using some standard thermodynamic relations 
(giving the partial derivatives of the pressure and of the rescaled chemical potential with 
respect to the temperature) the various kinetic coefficients can be determined, after some 
algebra: 


Bb{Pr,T) = T ttBipR), 

Au){p riT) = 2T‘^aB{p.R)i 


Arb = 


d 


TaB^Pi 


^b{Pr,T) = AArPb'V^ 


(2.19) 

( 2 . 20 ) 


where aRijlR) is an arbitrary function of the rescaled chemical potential. Note also that 
Ab is fully determined in terms of the coefficient of the anomaly and it is, in practice, only 
function of the chemical potential itself since, by definition, ErB = Pr- 

In summary, in a globally neutral plasma with an anomalous current, the relativistic 
second law implies that the non-anomalous current must contain magnetic and vortical 
contributions. If the plasma is not hypercharge neutral the form of the kinetic coefficients 
is subjected to a higher degree of arbitrariness since a second chemical potential must be 
introduced in the analysis (see the appendix of Ref. 0)- 
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3 High-frequency propagation 


3.1 Two-fluid AMHD equations 


The kinetic coefficients of Eqs. fl2.19p and fl 2 . 20 p can be redefined, for practical reasons, as: 


T) — SttT Cij(/i^,T) — 167rAi{/i^T, 

TaBipn) 


crb{^-r-,T) - 


dji 


R 


CRuj{I^Ri T) — 0. 


(3.1) 


The hypermagnetic and hyperelectric fields denoted, in what follows, by B and E obey the 
following set of equations 


V ■ E = 47rq'(?7,+ — n_), V x E = —dtB, (3.2) 

VxB = — n_vJ) + c^u — cbB + dtE , (3.3) 

where B is divergenceless (i.e. V-S = 0). The coefficients c^(J1r,T) and cb{JIb.,T) multiply, 
respectively, the vortical and the magnetic currents of Eq. (13.3p . The three-vector where u 
defines the total vorticity u = (p+a;+-l-p_a;_)/(p+-|-p_) and should not be confused with the 
frequency (denoted by hi in what follows). The energy densities of the charged species are 
denoted by p±. To establish a direct connection with the spectrum of conventional plasmas 
in the limit of vanishing kinetic coefficients we shall preferentially consider the situation 
where the charged species are massive, namely p± = n±{m± -|- 3T±/2) and p± = n±T± with 
T±/m± < 1 . 

The remaining coefficients cruj{jJb,)T) and CRBiER^T) of Eq. fl3.ip affect directly the 
evolution of the concentrations: 

dtn+ + V • (h+ v+) + -V ■ {cuj A+a;+) - -V ■ {cb X+B) = 0, (3.4) 

q q 

dfi- -I- V ■ (h_ V-) -V ■ {coj X-U-) + -V ■ {cb X-B) = 0, (3.5) 

q q 

difiR + V • {ur vr) + V ■ {crui cJ) — V ■ (crb B) = —AArE ■ B, (3.6) 


where X± = p±/{p+ + p_). Concerning Eqs. fl3.2p - fl3.3l) and Eqs. fl3.4l) - fl3.6p few comments 
are in order. If h+ 7 ^ h_ a second chemical potential py (corresponding to the hypercharge) 
can be introduced in Eq. fl2.10p . The global hypercharge neutrality of the plasma implies 
Py = 0. The peculiar velocities determining the currents obey the following set of equations: 

dtV- + (n_ ■ V)n_ = —q^^lE + x .B] -|- rc(n+ — vA — , (3.7) 

P- P+ 

dtv+ + (n+ ■ V)n+ = q^[E + v+ x B] + Vc^{v- - v+) - (3.8) 

P+ P+ P- 

dtVR + {vr ■ V)vr = 0, (3.9) 


where the pressure gradients shall be eventually neglected; Tc denotes the collision frequency 
determining the generalized conductivity in the single fluid limit. 
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3.2 Linearization of the two-fluid equations 


Equations fl3.2p - fl3.3p . fl3.4p - fl3.6p and fl3.7p - fl3.9p will now be linearized in the presence of 
the weak background magnetic held Bq with the aim of deriving the dispersion relations. 
The background held will be considered homogeneous: this means that the variation of Bq 
occurs over typical length-scale^ much larger than 1/cb- The huctuations of the various 
quantities will be introduced as follows: 

n±{t,x) = nQ + 6n±{t,x), nR{t,x) = ni + 6nR{t,x), B(t,x) = Bo + 6B{t,x), (3.10) 

while for the other quantities (i.e. v±{t,x) = 6v±{t,x), vjiit,x) = 6vR{t,x) and E{t,x) = 
6E(t,x)) the huctuations coincide with the held itself. In Eq. fl3.10p no and ni are, respec¬ 
tively, the uniform background charge and the uniform chiral concentration. The homoge¬ 
neous value of the chemical potential is related to rii and the kinetic coefficients will also 
be homogeneous. In the case of approximate thermal equilibrium the chemical potential 
can be related to the concentration as = HoUR/q where denotes the entropy density 
at equilibrium and where /xq is a numerical constant. Therefore if ub is perturbed around 
a homogeneous background the kinetic coefficients will also be, in the hrst approximation 
homogeneous. Thanks to Eq. (I3.10p the perturbed version of the evolution of the concen¬ 
trations can be written as: 

-I- uo(V ■ = 0, 5n'j^ -|- ?7,i(V ■ 5vr) = -AA-r. 5E ■ Bq, (3.11) 


where the prime denotes a derivation with respect to the time coordinate t. Since the kinetic 
coefficients are homogeneous in the hrst approximation, their contribution disappears from 
Eq. fl3.1ip . With the same notations Eqs. fl3.7|) . fl3.8p and fl3.9p imply instead: 


= ± 




m± . 


6E -|- 6v± X Bi 


SvL = 0, 


(3.12) 


where T^ has been neglected but it will become relevant at low frequencies, as we shall see 
later. Finally, after inserting Eq. fl3.10p into Eqs. fl3.2p and fl3.3p we obtain: 


V-6E = 47iq{6n+-6n_), V ■ 55 = 0, 
V X 6B = 6E'+ Att qno{6v^ — Sv_) 


V x6E = -5B', 


- cbSB + 


A+V X -|- A_V X 


(3.13) 

(3.14) 


From Eqs. 03.121) the equation obeyed by 6u± can also be deduced and they are = 
Eq[—6B' -|- V X {6v± x Bo)]/m±. Recalling the standard vector identitie^ the equation for 

® In section |4] we shall specifically discuss also the opposite limit where Bq varies appreciably over typical 
lengths L < 1/cb and we shall see that, in this case, the background solution belongs to the class of Beltrami 
fields. 

^In particular we recall that V x (a x 6) = [a(V • b) — 6(V ■ a) + (b ■ V)a — (a • V)&]. 





































5oj± can also be expressed as: 


5oj'^ = ± 


q. 


m± . 


-6B' - Bo{V ■ 6v±) + {Bo ■ V)6v± 


(3.15) 


From Eqs. fl3.12p and fl3.14|) the relevant dispersion relations and the associated refraction 
indices can be obtained by treating separately the motions parallel and perpendicnlar to the 
magnetic held direction. 


3.3 Appleton-Hartree determinant 

While in conventional plasmas the Appleton-Hartree dispersion relation has been extensively 
discnssed in the literatnre mm, the AMHD eqnations linearized in the two-hnid limit 
contain vortical and magnetic cnrrents. The Laplace transform of Eq. fl3.14p implies the 
following eqnation 

(V X SB)q = -iQssi^) ■ 6 Eq - cbSBq + ic^V x • SEq], (3.16) 


where H is the freqnency (not to be confnsed with the total vorticity). In Eq. fl3.16p ^^(H) 
and £^(H) denote, respectively, the standard and the vortical components of the dielectric 
tensor. The explicit form of ^^(H) and £u(H) can be fonnd in appendix taking then the 
cnrl of Eq. fl3.13p and nsing Eq. fl3.16p we obtain the following eqnation; 


V X (V X 6En) = nh,{n) • 6En - csV x H V x • 6En]. (3.17) 


We can now go to Fonrier space and write Eq. fl3.17p as: 


— k X k X SEj^^ = £s(H) • SEj:^ — icsk x SEj:^ — ic^^Qk x [£„(H) • SEj:^]. (3.18) 

We can therefore introdnce the refractive index[^ n satisfying n = k/Vt where k = |fc|; 
choosing the coordinate system as fc = (0, nVL sin 0, nVL cos 9) we can obtain from Eq. fl3.18p 
the following Appleton-Hartree matrix: 


/ [l-^ + c^^c(0)] 

+ ^c(0) -I- ^eoc{e)] 


+ ^cW + ^^3c(0)] 

-s{e)c{e) + ^£ 4 s( 0 ) 


-s{9)c{d) . (3.19) 

[,2(0)_i^] ) 


The above matrix rednces to the standard form of the Appleton-Hartree matrix in the limit 
Ctj —)■ 0 and cb —t 0. 

®The refractive index cannot be confused with the concentrations denoted by n± and ub] their homoge¬ 
neous values no and ni carry specific subscripts so that the notations are clearly established. 
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4 Dispersion relations 

The determinant of the Appleton-Hartree matrix obtained in Eq. fl3.19p leads to the following 
expression: 

sin^ 6 '(e|| - n'^) [71^(81 + Sr) -28^81] - 2 cos‘^ ee\\{n^ -£L)(n^ - Sr) 

+2n®[c|/B(£, n, n, 9) + CRgBis, n, 9) + clf^{e, hi, n, 9) + ff, n, 9) 

c| hi{e, 12, n, 9) + crc^ h 2 {e, 12, n, 0) + cl cr h{e, 12, n, 0)] = 0. (4.1) 


Equation fld.ip is written in terms of the 7 functions explicitly reported in Eq. flA.711 of 
appendix These functions have a specihc dependence upon the dielectric tensors; with a 
collective notation such a dependence has been indicated by e. The notations followed in 
Eq. fld.ip imply that c% multiplies fs, cl multiplies fuj] gs and g^ multiply, respectively, cr 
and c^; the three functions hi, h 2 and multiply instead the mixed products. Finally both 
in Eqs. fl4.ip and in Eq. flA.7p we have introduced Sr = {si+ 62 ) and Sr = (^i — 62 ) given 
by: 


£ 1 ,( 12 ) = 1 


12 ^ 


(12 + 12 b _)(12 — 12 b _|_)’ 


^i?( 12 ) — 1 


12 ^ 


(12 + 12 ^_ i _)(12 — 125 _)’ 


(4.2) 


where 12^ = (12p_,_ + 12p_). When cs = = 0 the magnetic and the vortical currents are 

absent from the two-fluid AMHD equations and Eq. (14.ip implies the standard result [6]: 


sin^0 




1 

£11 


1/1 


+ 


.2 Ker Sr 




= cos^9[ ---- 


£l 








(4.3) 


The dispersion relations for a wave propagating parallel (i.e. 0 = 0) and perpendicular (i.e. 
9 = 7r/2) to the magnetic held direction can be obtained easily derived from Eq. (14.3p . If 
0 = 0 Eq. (14.31) reduces to (n^ — £r)(u^ ~ £l) =0 while for 0 = 7r/2 Eq. (14. 3 p implies 
{in? — £m)[?7,^(£l + £r) — 2£l£r] = 0. These dispersion relations give therefore the conventional 
result^ which will be generalized hereunder. 


4.1 Free-field propagation 

In the absence of magnetic held there is no preferred direction and the dispersion relations 
follow from Eqs. (14. 3 p and (lA.7p by setting all the Larmor frequencies to zero. In this case 
eR = Sr = e\\ and the dispersion relations stem from the following two conditions, namely: 

£||(f2) = 0, {r? — £\\)VL^ ncR = (4.4) 

®Along 0 = 0 we thus obtain usual dispersion relations for the two circular polarizations of the electro¬ 
magnetic wave, i.e. r? = er and r? = el, while along 0 = 7r/2 we have the dispersion relations for the 
“ordinary” (i.e. r? = Ey) and “extraordinary” (i.e. r? = 2er,El/(er -I-El)) plasma waves. 
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Equation fl4.4p demonstrates that the vortical current does not contribute to the dispersion 
relations in the free-held case: is absent from Eq. fl4.4l) since the two-huid effects cancel 

in the total vorticity. This cancellation is either exact (as in the case of the free-held propa¬ 
gation) or approximate (as we shall see later in the presence of the magnetic held). Indeed, 
as it can be explicitly verihed from Eqs. flA.ip . flA.Sp and flA.bp . et,(r2) —)■ 0 when Bq —)■ 0: 
in the limit i?o —>■ 0 the vorticity of positively and negatively charged species is balanced so 
that the net total vorticity vanishes. 

The dispersion relation £||(hl) = 0 implies hi = hip. This wave does not propagate since 
its group velocity vanishes and these are nothing but the electrostatic plasma oscillations 
[0]. The solution of the second equation in Eq. fl4.4p is insteaci^ 


” = =^2S + V^" 


+ 




9? A9?' 


(4.5) 


Equation fl4.5l) implies also + k'^ ^ kcB] these modes are propagating but only 

ahected by the magnetic current, as previously remarked. The birefringent nature of the 
dispersion relations will be discussed a bit later since this free-held ehect may interfere with 
the presence of the background magnetic held. 

If Cs —?■ 0 and —)■ 0 we have that n ^ 0 whenever one of the following three possibility 
are separately verihed e:||(hl) = 0 or £1,(14) = 0 or even £h( 14) = 0. The frequencies arising 
from the previous conditions are cut-ohs because, for given equilibrium conditions, they 
dehne frequencies above or below which the wave ceases to propagate at any angle (fc —)■ 0 
for hnite 14, i.e. Up = 14/fc —)■ cxd). This is what happens, in particular, with the dispersion 
relation of Eq. fl4.5p . Let us hnally remark that the remaining two cut-ohs stemming from 
the conditions £l( 14) = 0 and £_r( 14) = 0 in Eqs. fl4.2|) are given, respectivey, by: 


14r — Y^14p + (14s+ + 14s-)^/4 — (14^+ -I- 145_)/2, (4.6) 

14l = Y^14p -l- (14b+ + 14_b-)^/4 -I- (14ij+ -I- 9b-)/2. (4.7) 


4.2 Parallel propagation 


Taking the limit 6 ^ —)■ 0 in Eq. fl4.ip and recalling the results of Eq. flA.7p we obtain: 

£\\{ncB + [n^ + nc^ (ss + £4 ) - £k]14}{ucs - [n^ + nc^ (-S 3 + £4 ) - £l]14} = 0. (4.8) 


If £||(14) = 0 we go back to the case of electrostatic oscillations. Therefore, assuming £||(14) 7 ^ 
0, Eq. fl4.8p implies that the standard dispersion relations are modihed as: 


n 

n 


1 

1 


— Cb — Ccj(S 3 -h £ 4)14 ± 4£jj14^ -|- [cb + C^14(£3 -1- £ 4 )]^ 5 
Cfl + C„(Es - Ei)n ± + |cb + C„S1(£3 - £ 4 )]^ . 


(4.9) 

(4.10) 


^°The positive square root has been chosen in Eq. (14.5p in order to get > 0; we consider only positive 
VL since solutions with 17 < 0 simply correspond to waves travelling in the opposite direction. 
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Thus the dispersion relations for the generalized L-mode and -R-mode are given, respectively, 
by: 


+ k[cB + Cuj{s3 + £ 4 )^^], 
e:L(hl) = P - k[cB + - £ 4 )^!]. 


(4.11) 

(4.12) 


In the high-frequency limit (i.e. formally hi —?• cxo) we have that C(^(e :3 ± £ 4)12 —)■ 0 since, 
from Eqs. (IA.5ll - flA.6ll . we have: 


Ctj(£3 + £ 4 )^^ — 
Clj(£3 ~ £4)^^ = 


qc^ 


(m+ -|- m_) 
qcuj 

(m+ -|- m_) 


Q 


Qb--^} 

n 


+ 


Q 


+ + fi¬ 
ll 


fl — 11B + H T H-B —- 


(4.13) 

(4.14) 


The results of Eqs. fl4.11l) - fl4.12l) and fl4.13p - fl4.14p demonstrate, once more, that in the 
high-frequency limit of the spectrum the magnetic current dominates against the vortical 
current. For intermediate frequencies (i.e. as soon as we reduce H) the terms containing the 
natural frequencies of the plasma come then into play so that for the R and L modes the 
corresponding dispersion relations become: 


H — k kcB T 


H'll 


= k"^ — kcB + 


(H-Hb-)’ 

n^n 

(n -|- Hb _) ’ 


R — mode, 
L — mode. 


(4.16) 


As in the standard case, the phase velocity of the .R-mode is greater than that of the L-mode. 
In Eq. fl4.15p we assumed m_|_ > m_ and therefore Hb-i- < Hb-- In fhe limit A; —)■ 0 the 
R-mode cut-off occurs above Up while the L-mode cut-off occurs below Up (i.e., recalling 
Eqs. (14. 6 p and (14.71) . H —)■ and H —)■ Hl). In the low-frequency limit Sr and El coincide 
to leading order in (11/11 b_|_) and in (11/11 b_) since 


112 

Ihn ^(H) = lim £^(11) ^ 1 + -^— 

El—^O Q—>-0 SiB-f-ilB — 


= 1 + 


Bn 


2 ’ 


Va = 


47rno(m_|_ -|- m_) 


(4.16) 


where va denotes the Alfven velocity of the system. In the low-frequency limit the dispersion 
relations for the .R-mode and for the L-mode are, respectively. 


H" = 




11^ = 


1 + ni 

Va 

1 + vl 

2 //'I I „,2 ■ 


k^ + k 


cb + 


qcc. 


n 


m V11b + 


+ 


H 


k^-k 


Cb — 


qcu. 


H 


+ 


Hb-/J 
H 


m VHb-i- h 


B-/ J 


(4.17) 

(4.18) 


since -C 1 we also have Va/{)- + Va) — va 
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Having determined the dispersion relations in the case of parallel propagation, the Fara¬ 
day rotation rate can be easily determined with the standard procedure. The generalized 
Faraday rotation angle experienced by the linearly polarized radiation travelling parallel to 
the magnetic field direction can be obtained as 


Ad) = 


_ n 
“ 2 I 


+ CcjEs + 1/ Sl -|- 


^ + S4) 


^LO / \ 

^ + y(£s-£4) 


— \ £r + 


AL, 
(4.19) 

where AL is the distance travelled by the signal in the direction parallel to the magnetic 
field direction. It is interesting to compare the contribution of the terms depending upon cb 
and those depending upon the background magnetic field intensity, i.e. the terms appearing 
in the squared brackets. Recalling the expressions of {er, El) we have that flB+ ^ 
and Qp- (always assuming 3> m_). In this case Ad)/AL interpolates between 

the standard result {Qb-/‘^){^p-/^)'^ (valid when cb —t 0) and the constant rotation rate 
cbI‘2 (valid when Rq 0 as in the case of free-held propagation). As it can be explicitly 
verihed the is subdominant at high frequencies and can be neglected. 


4.3 Orthogonal propagation 


By setting 6* —>■ 7r/2 in Eq. (14. Ih and recalling the results of Eq. 
simple equation: 


Tj) we obtain the following 


n^{EL + — n^{c%{EL + Er) + CbCi^[E4{—El + Er) + Ez{El + 

+ [^||£k + ^l(^|| + 2e_R)]h2^} -I- 2elErVl^e\\ = 0. (4.20) 


The solution of Eq. (14.201) can be obtained by hrst solving in terms of n^. The result is 


+ 

M{E,n) = 
+ 


j{e, n) ± ^m{e, n) 

2{er + El) 

C^b{^L + £r) + Cb Ci_j [E/^^—El + Er) 4- E^{El + £_R)]f^ 

+ Cl{e\\ + 2£ij)]r2^, 

— ‘&C\\CrEl{cL + -|- {c\{El + Er) 

cbCu][eii{er — El) + + £R)]f2 -|- [e||£_R + £l{,c\\ + ‘2er)]VL^)'^. 


(4.21) 

(4.22) 

(4.23) 


Equation (I4.2ip in the limit —)■ 0 and cs —>■ 0 reduces to the ordinary mode if we choose 
the plus (i.e. in? = ey) and to the extraordinary mode (i.e. r? = 2 erElI{eb + £l)) if we 
choose the minus. In the high-frequency limit the terms multiplying the vortical current are 
always negligible as already remarked in the case of the parallel propagation. The phenomena 
related to the oblique propagation will not be specihcally discussed. 
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4.4 Spectrum around a hypermagnetic knot 

Introducing the three mutually orthogonal unit vectors a{z,p), b{z,p) and z dehned in ap¬ 
pendix [HI we can consider the modes of fluctuation of the hypermagnetic held around a fully 
inhomogeneous background BQ(t,x), namely: 

B(t, x) = Bo{t, x) + 5B(t, x). (4.24) 

Since the background solution is not uniform we can align Bo along a{z,p) and write that 
Bo{z) = Bq a{z, p). The background equations are solved by setting p = —cb (since V xBq = 
pBo). As in the homogeneous case the velocites vanish on the background solution. For 
L < 1/cb the background held is homogeneous and the previous analyses apply. For typical 
length-scales larger than the scale of spatial variation of hypermagnetic knot (i.e. L ^ l/c^) 
there are two separate possibilities for the perturbed velocity held: either (5F || Bq or 6 v B Bq. 
These two cases will now be separately examined. 

The case of parallel propagation mirrors exactly the one already discussed in the case of 
uniform held. If we assume that 5F || Bq the dispersion relations follow from 


{[/c^ — f2^e||(r2)]^ — c%k‘^}e\\{Vt) = 0. 


(4.25) 


The parallel dielectric tensor is £||(f2) = 1 — np/[f2(f2 -|- iFc)] where the correction coming 
from the collision rate has been added for immediate convenience. The solution e:||(f2) = 0 
gives, as before, the electrostatic wave. The solution of = 0 gives, 

respectively, a high-frequency and a low-frequency branch. The high-frequency branch has 
the same dispersion relation of the free-held case, namely 12^ ~ 12^ -|- =F kcB- The low- 
frequency branch is instead derived from the explicit form of the dispersion relation written 
as: ^ 

^ = ^(12 + *F,) - M ± ^(12 + *F,); (4.26) 

neglecting the hrst term at the right-hand side of the previous equation (which is unimportant 
at low frequencies) we have that: 

^ _ zFc(l =F CB/k) 

1 + ni/k‘^TCB/k' 

The low-frequency mode, in which the conducting current dominates over the displacement 
current, has no counterpart in vacuum. In the low-frequency mode, a small electric held 
proportional to F^ exist to give the necessary current parallel to the magnetic held. In 
the limit F^ —)■ 0 both the electric held and resistivity vanish and the low-frequency mode 
becomes the force-free held. As expected the same phenomenon occurs in the absence of 
magnetic and vortical currents [2011^ . 

In the case of orthogonal propagation the huctuations of the hypermagnetic held and of 
the velocity can be decomposed by using the gyrotropic basis of appendix jB) 

SB{t,z) = 6 Bi{t)b{z,p) + 6 B 2 {t)z, 6 'if^\t,z) = 6 vi^\t)b{z,p) + 6 v 2 ^\t)z. (4.28) 
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For a generic velocity fluctuation orthogonal to Bq we have SvxBq = [Bo{6v-z)b—Bo{6v-b)z]; 
the solutions for 6 v['^\t) and 5^2^^ (t) can then be expressed as: 


54^^^) = 


q 


The dispersion relations in this case are given by 

Cb Cuj 


6Ei + D.b ± bE2 
±in6E2-nB±SEi 




2Q 


^l(^3 ~ ^ 4 ) + £r{£3 + £ 4 ) 


= 0 . 


(4.29) 


(4.30) 


In the high-frequency limits dehned by Eqs. fl4.13p - fl4.14p . Eq. fl4.30|) is satished if EbSl = 0 
which is verihed when either el or er are vanishing. Equation fl4.30p leads to vanishing group 
velocity in the high-frequency regime: the corresponding modes are then not propagating. 
The proper frequencies dehned by these equations have been already derived in Eqs. fl4.11l) 
and fl4.12p . 


5 Single fluid description and its implications 


The two-huid equations can now be combined with the purpose of deriving the effective single 
huid description valid for sufficiently large length-scales and for frequencies much smaller 
than flp and ^Ib±- The one-huid variables are the total current J = q (n+, — n_v_), the 

bulk velocity of the plasma v = (m+n+ -|- rri-V-)/{m+ + m_) and the total mass density 
Pm = (m+n+ -|- m-n_). In the globally neutral case J and pm become, respectively, J = 
no(hV — vJ) and pm = nQ{m^ + m_). Summing-up Eq. fl3.7l) (multiplied by m+) and Eq. 
fl3.8p (multiplied by m_) the evolution equation for the bulk velocity of the plasma is 


Pv 


dtv + V ■ VF 


= J X B — VP -I- 


(5.1) 


where the shear viscosity contribution, labeled by r^, has been added for convenienci^. 
Taking the difference of Eqs. fl3.7l) and fl3.8p (and assuming m+ > m_) the generalized 
Ohm’s law can be written as: 


dJ+rJzE 


dvr 


(^E + vxB + 


Vp 

quo 


Jx B 
quo 


(5.2) 


where global neutrality has been assumed. Note that we have also kept the thermoelectric 
term (depending on the pressure gradient of the lightest charge carriers) and the Hall term. 
Since we shall mainly consider the case of homogeneous pressures the thermoelectric term 
will be neglected; the Hall term is a higher order contribution, as we shall argue. 

the total pressure does not vanish Eq. dSH) is modified as follows ft [n; F| -I- (n • V)[r(;'(;] -I- vV ■ [wn] = 
—VP-I-J X B + rjV^v where w, as already discussed, is the enthalpy density. 
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In the globally neutral case the single fluid equations stipulate that E, B and J are all 
solenoidal (i.e. V--E' = V- -B = V- J = 0). A fourth possible solenoidal vector is the bulk 
velocity of the plasma v. Indeed, since the evolution of pm and pq = g(n+ — n_) is given b}@ 

dtpm + V ■ {pmv) = 0, dtpq + V ' J = 0, (5.3) 

the incompressible closure V • n = 0 will be adopted; consistently with the incompressible 
closure pm will be considered homogeneous, at least in the hrst part of this section. A 
full discussion of other possible closure (such as the ones conventionally adopted in conven¬ 
tional plasmas) is desirable but beyond the scoped of this analysis. The remaining one-fluid 
equation containing the vortical and the magnetic currents, can be written as: 

V X 5 — dtE = 4:71 J + — cbB, uj = V x v. (5.4) 


Since the one-fluid description involves the lowest branch of the spectrum we can neglect 
the displacement current that becomes relevant only for the electromagnetic propagation. 
For the same reason we can neglect the time derivative in Eq. fl5.2|l . i.e. dtJ -C ^pE. 
Consequently Eqs. 115.211 and 115.4p in the low-frequency branch of the spectrum become 

V X B = AttJ + c^cJ — CbB, E = J j a — v x B. (5.5) 


Recalling that V xE = —dtB, Eq. 115.511 can be used to obtain an equation that is reminiscent 
of the magnetic diffusivity equation, namely 


-> -> -> 1 -> 1 
dtB = 'V X {v X B) + - h -—V X — - —V x {cbB). 


47ia Alia 


dvrcT 


(5.6) 


Introducing the vorticity uj into Eq. 115.111 and dividing both sides of the equation by p, 
we obtain 

P ^ 


o J xB ^ 

OtV + u X V = -V 


Pri 


-pr 


H” ^kin^ 


(5.7) 


Taking the curl of Eq. 115.7p the evolution equation of the vorticity becomes: 


dtu = 'Vx{vxu}) + 


V X (J X P) 2- 

-1" Xl. 

pm 


(5.8) 


The most interesting solutions of the one-fluid equations will involve the situations where 
the vortical and the magnetic currents play the dominant role. However, before turning the 
attention on these classes solutions it is useful to remark that the equilibrium solutions of 
the plasma at rest (i.e. n = 0) are simply given by 


VP = J X B, AttJ = I V X P — c^uj + cbB ). 


(5.9) 


Note that the global neutrality implies Pq = 0 and Eq. (15.311 demands V • J = 0 in full agreement with 
the solenoidal nature of the total Ohmic current. 
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From Eq. fl5.9p it is immediate to obtain the following three identities 

AttVP = {V X B) X B, (S-V)P = 0, (J ■ V)P = 0. (5.10) 

The relations of Eq. fl5.1Up are not explicitly modified by the presence of vortical and 
magnetic cnrrents. The two last conditions in Eq. fIS.lOp dehne the so-called magnetic 
snrfaces: the pressnre gradient vanishes along the lines of magnetic force and along the 
current lines. We conclude that neither the vortical not the magnetic current affect directly 
the equilibrium solutions. 


5.1 Bulk velocity parallel to the magnetic field direction 

The single fluid AMHD equations admit various solutions that have no counterpart in the 
case of ordinary MHD. Consider first the situation where the hypermagnetic magnetic field 
and the velocity are parallel and have non-vanishing magnetic and kinetic gyrotropy, i.e. 

V X B = 0, V ■ V X V = Pv{t)v‘^, B -V X B = pB{t)B‘^. (5.11) 


The simplest way to realize the situation described by Eq. fl5.11|) is to require that v 
and B are both Beltrami-like helds (see appendix [B] for this terminology) characterized 
by V X n = Pv{t)v and V x P = pB{t)B. Moreover, since n x P = 0, it is natural to require 
that Pv{t) = ps(t). From Eq. fl5.5p the total current J can be easily determined; the Ohmic 
electric field is then given by: 


tj — --- B --- OJ. 

dvrcr dvrcr 


(5.12) 


From Eq. fl5.6p the hypermagnetic held is obtained by solving the following equation 


dtB 


VBit)[VB{t) + CBit)] p 

Anait) 


Cujjt) 

A7ra(t) 


PB(t)S, 


(5.13) 


where u}{t,z) is the solution of Eq. fl5.8l) . Thanks to the symmetries of the problem the 
solution of this equation is given by uj{t,z) = a;o(-s)exp [— Jq PB{'t')^kin{'t')dt'], where uJo{z) 
is the initial vorticity which can also be written as uo{z) = pb{0)vo{z). Equation 05.131) can 
then be solved in general terms. However, recalling that cb( 0 and Ccj{t) are explicit functions 
of time but they depend on the rescaled chemical potentiate since p_b(^) is arbitrary we can 
choose ps(t) = —CBit). In this case, the solution of Eq. 05.13P shares the same properties of 
the general solution but it is mathematically simpler: 


^(*'")="»(")r*'^i') 


(0 -/■' 

-—e Jo 




(5.14) 


this discussion we shall keep the time-dependence in the kinetic coefficients even if, strictly speaking, 
cb and c^, may depend on the temperature and of the chemical potential but they are constant in time. 
However, in curved backgrounds a mild breaking of conformal invariance may induce a time dependence 
which is, however, not central to the present analysis. 
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The results of Eq. fl5.14|) describe the generation of the hyper magnetic held thanks to some 
initial vortical current. To deepen this question let us assume that c^j and cb are both 
constant. Equation fl5.14p can then be solved and the result is 


B{z,t) 


hb(^) 

dvr l^kin ^ 




(5.15) 


where we used that uo{z) = —cbVo{z) when cb is constant in time and cb(0) = cb- This result 
is also valid for a relativistic equation of state (i.e. w = 4p/3) provided the incompressible 
closure is consistently adopted and can be easily generalized to curved backgrounds. While 
these generalizations are not germane to our theme it is worth to emphasize that in the 
limit t ^ oo the suppression of the magnetic held is controlled by Airavkin which is nothing 
but the Prandtl number given as the ratio of the magnetic and of the kinetic Reynolds 
number [2]. The Prandtl number is roughly independent on the temperature. For instance 
at the electroweak epoch m we would have that Ukin — while a ~ T/a' where 

a' = /{An). Recalling the results of Eqs. fl2.19p and fl2.20p we therefore have that 


lim B{t, z) 

t^oo ^ ' 


C^{T) 

An 




c^(T) = 2T^a(/r^), 


(5.16) 


where aBiJAn) is the usual arbitrary function of the rescaled chemical potential = fin/T 
(see Eqs. (I2.19p and (I2.20p ). The suppression due to the conductivity is therefore eliminated 
and what is left is a milder suppression 0{a'^). In the same limit the hypermagnetic current 
turns out to be more suppressed than the vortical current. As long as the inverse of the 
Prandtl number scales as a'^ the previous discussion is generally valid and this is what 
happens in the case of the electroweakplasma [lO] (see also [H] for specific estimates of the 
conductivity in the electroweak phased . 


5.2 Bulk velocity orthogonal to the magnetic field direction 

If the hypermagnetic field and the bulk velocity of the plasma are orthogonal (i.e. v- B = 0), 
employing the helical basis of appendix [B] the hypermagnetic held and the velocity held can 
be written as: 

B{t,z) = Bo{t)z + Bi{t)a{z,p), v{z,t) = v{t)b{z,p). (5.17) 

From Eq. fl5.5p the AMHD current becomes: 

= ^^[cB+p]a{z,p) - ^pv{t)b{z,p) + (5.18) 

By analyzing the structure of the evolution equation of the magnetic held and of the vorticity 
it emerges that the system is consistent provided dtBo = 0 (i.e. constant magnetic held along 

^^Similar kinds of considerations s can also be developed in the case of a strongly interacting plasma as 
long as the same scaling occurs m- 


18 





















z) and provided = 0. In this case the conpled evolntion of the vorticity and of the magnetic 
field obeys 


duj 

dt 


PCrBo 

1 - 


[cb + p]Bi, 

4:TTpm 


dBi 

dt 


uBo - 


pBi 

dvrcr 


[cb+p1, 


(5.19) 


where u{t) = pv{t) and u{z,t) = u{t)b{z,p). The eqnations can be diagonalized with a 
specihc choice of the coordinate system. The simplest and most convenient one is p = —cb] 
in this case the two eqnations can be combined by differentiating once Eq. fl5.19l) . The resnlt 

is 


d‘^Bi 

dB 


+ c%v\Bi 


0 , 


Bo 

^/Wp, 


(5.20) 


where va is the Alfven velocity. This solntion has been swiftly presented in Ref. [7] and 
recently rediscovered in [23]. Eqnation fl5.20p describes the AMHD analog of the non-linear 
Alfven wave. The anomalons Alfven wave has been already discnssed in section [3] as a 
low-freqnency limit of the two-flnid eqnations. 


5.3 Fully non-linear evolution and baryogenesis 

So far we considered small flnctnations of the chiral concentration aronnd an otherwise 
homogeneous value denoted by rii in section [21 In the opposite case the AMHD equations 
imply a specihc relation between the concentration (or the chemical potential) and the 
topological properties of the hypermagnetic helds. To illustrate this point we will show 
that close to equilibrium the chemical potential is determined not only by the magnetic 
gyrotropy but also by the total vorticity of the plasma. Hypermagnetic held conhgurations 
with non-vanishing gyrotropy have been used to model the generation of the baryon or 
lepton asymmetry [TO] (see also [251126112Z]). Consider, therefore, the evolution equation of 
the chemical potential which can be written as 

dt-pn + = -^ArE ■ R, (5.21) 

where T is the perturbative rate of he chirality hip processes (in the case of [10] it is deter¬ 
mined by the scattering of right electrons with the Higgs and gauge bosons and with the 
top quarks because of their large Yukawa coupling). In Eq. (I5.2ip we also used the follow¬ 
ing general relation relation = PqUr/^, where po is a numerical factor depending on the 
specihc features of the plasma while ? is, as usual, the entropy density. 

To compute p^ in the proximity of an equilibrium situation we need to deduce the 
hyperelectric held. Recalling then Eq. f|5.5p . the hyperelectric held ca be related to the total 
Ohmic current so that Eq. (15.211) will become 

dtPB + T/Ij; - ^ (^b{Pr) ^ ^ ^ iq PoTPr ^2 _ V X B. (5.22) 

cr<^ aq naq 
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We shall now choose aBiJIji) = if aB{pji) ^ J^b fh® evolution equation of the concen¬ 
tration we should add a further term proportional to V • [cb,b{'Pr)B]. This term vanishes in 
the case aBipR) oc JIr: crb{Pr) contains the derivative of qr^JIr) with respect to JIr and it is 
therefore constant. In more general situations Eq. fl5.22p will just contain a supplementary 
contribution of the type B ■ VJIr. 

The rescaled chemical potential enters the inhnitely conducting limit (see appendix [B]) 
and the smallness of the particle asymmetries is the rationale for the minuteness of the 
rescaled chemical potentials in approximate thermal equilibrium. At equilibrium, we can 
determine JIr from Eq. (15.2211 and the result is given by: 


/B ■ V X B 

1777; [T + Ts - r^] ’ 


(5.23) 


While r is the perturbative chirality flip rate, the other terms can be understood as rates 
stemming from the hypermagnetic current and from the vortical current and they are 


Tb = 

^cr 


^ _ 8T7o ^ ^ B 
—- Aru ■ B. 

TT a<^ 


(5.24) 


In the case of right electrons (see |T0]) Ar = —5f'^?/|./(647r^) where g' denotes the gauge 
coupling and ur = —2 is the hypercharge assigment of the right electrons. In the same 
situation we have JIr = fioUR/g and fio = 877r^ iVpy f/220, where W// is the effective number 
of relativistic degrees of freedom of the systeno. If the plasma is hypercharge neutral the 
value of the chemical potential can be estimated from the asymmetry in the case where all 
the standard model charges are in complete thermal equilibrium. If all the asymmetry is 
attributed to the right electrons (which is, in some sense, the most favourable situation) 
then Jlj^ = /220) Nesinn/^^) where We = 106.75. With these specihcations Eq. (I5.23p 

becomes 

783 a' B-VxB ^ 783 a'^ B^ 

~ ~88 7raT3 [r + Ts - Tj ’ ^ “ 22 7rVT3’ ^ 

where a' = g'‘^/{ATr). Equation fl5.25p coincides with the previous results (see e.g. Eq. (6.15) 
of the last paper of [10] and see also [25l [26l |2^) in the limit T^ —)■ 0. The results of Eq. 
fl5.25p show that the hnal value of the chemical potential depends on the properties of the 
flow entering the dehnition of T^. In summary we can say that the magnetic currents and 
the vortical currents can affect a number of processes such as the formation of the baryon 
asymmetry or the dynamics of the electroweak phase transition. Similar kinds of effects can 
be expected in the case of strongly interacting plasmas where the magnetic gyrotropy can 
also determine the properties of the chemical potential. 

^^There have been a number of suggestions for possible roles that the abelian hypermagnetic Chern-Simons 
term might play in cosmology. One of them is related to the observation that right-handed electrons, which 
do not take part in weak interactions and also have a very small Yukawa coupling, are practically decoupled 
from the thermal ensemble above temperatures of about 10 TeV. 
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6 Concluding remarks 


The dispersion relations of anomalous magnetohydrodynamics are affected by the vortical 
and the hypermagnetic currents. The vortical currents do not impact on the high-frequency 
branch of the spectrum but the opposite is true at lower frequencies where new solutions 
describe the simultaneous presence of hypermagnetic knots and fluid vortices. These parity- 
odd configurations carry, respectively, hypermagnetic and kinetic gyrotropy. The physical 
properties of the system roughly interpolate between the features of conventional chiral 
liquids and the results valid for cold electromagnetic plasmas. While chiral currents are 
anomalous and do not contribute to entropy production, vector currents are associated with 
the generalized Joule heating. 

When chiral and Ohmic currents are simultaneously present the second law of thermody¬ 
namics constrains the kinetic coefficients. The hypermagnetic, vortical and Ohmic currents 
affect the evolution of the gauge helds and determine the hyperelectric held of the plasma. 
In anomalous magnetohydrodynamics the perfectly conducting limit is well posed and the 
hypermagnetic helicity of the knots is strictly conserved, as it happens in the case of conven¬ 
tional plasmas. The hypermagnetic currents are then completely washed out in the perfectly 
conducting limit and strongly suppressed when the conductivity is large but hnite. Close to 
thermal equilibrium the concentration of the chiral species and the corresponding chemical 
potential will depend not only on the hypermagnetic gyrotropy but also and on the vortical 
currents. 

In summary the evolution equations of anomalous magnetohydrodynamics offer a minimal 
theoretical framework where the interplay between conduction currents and chiral currents 
can be quantitatively analyzed. It is therefore fair to say that the results derived here 
complement and extend some of the present and earlier strategies aimed at an improved 
understanding of chiral liquids when generalized Ohmic effects cannot be neglected in the 
evolution of hypermagnetic and hyperelectric fields at hnite fermionic density. 
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A Generalized Appleton-Hartree equation 


A.l Explicit form of Eg and 


We are going to give, in what follows the explicit form of the dielectric tensors appearing in 
section 131 The matrix form of and is given by: 


(A.l) 



/ £i(ll) 

i£ 2 (ll) 



! £ 3 ( 12 ) 

— i£ 4 ( 12 ) 0 \ 

£s(ll) = 

-i£ 2 ( 12 ) 

£i(ll) 

0 , 

£^( 12 ) = 

i£ 4 ( 12 ) 

£ 3 ( 12 ) 0 j 


V 0 

0 

£||( 11 ) ) 


V 0 

0 0 / 


02 


^2 

P- 






B- 


“p- 


‘‘B- 


Q 


B- 


where £ 2 ( 11)5 £ 3 ( 11 ), £ 4 ( 11 ) and £||(11) are dehned as: 

£l(fl) = 1 

^ll(ll) = 1 

^ 3 ( 11 ) = 

£ 4 ( 11 ) = 


“P- 






02 
^‘p 4- 

qVt 


^^p - 
fl2 ’ 


(m+ + m_) 

q 

(m+ + m_) 


-f]2 


-f]2 


Q 


B- 




+ 


VL 


B- 


-ff2 


(A. 2 ) 

(A.3) 

(A.4) 

(A.5) 

(A. 6 ) 


Both £ 3 ( 12 ) and £ 4 ( 12 ) have dimensions of an inverse frequency squared; £i(12), £ 2 ( 12 ) and 
£||(12) are instead dimensionless. The frequencies appearing in Eqs. flA.2D - flA.4D and flA.Sp - 
(lAl) are the plasma and the Larmor frequencies associated with the charge carriers of 
both signs, i.e. f2p± = y/A'iiq^no/m± and 12 b± = qBo/m±. To compare the dispersion 
relations with the standard situation of cold plasmas we must bear in mind that the ratios 
of the plasma and Larmor frequencies are related to the inverse ratio of the masses, i.e. 
12p+/f2p_ = VIb+/^b- = m_/m+. 


A.2 The seven function 


The generalized form of the Appleton-Hartree equation (see Eq. (14.11) ) depends on 7 functions 
whose explicit form is given by: 

/s(£, 12, n, 0) = cos^ 9 + 


n^f22 




gB{e,n,n,e) = cos 0, 

POS^ f) ( 

fu}{E, 12, n, 9) = —(£3 — £ 4)(£3 + £ 4 )£|| + 
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sin 6*£3 — (£3 — 2e2) sin 2 0 |, 
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g^{e, n, n, 9) 


hi{e, f2, n, 6) 
h 2 {e, Q, n, 6) 
hsie, n, 9) 

The functions reported in Eq. (IA.7P determine, through Eq. fl4.ip . the form of the dispersion 
relations when the hypermagnetic and the vortical currents are simultaneously present in 
the anomalous magnetohydrodynamics equations. 

B Hypermagnetic knots and Beltrami fields 

In the resistive approximation, the hyperelectric and the hypermagnetic helds are not exactly 
orthogonal and the nature of this misalignment is crucial both for the generation of the 
baryon asymmetry and for the chiral magnetic effect. In AMHD the induced hyperelectric 
held stems directly from the approximate form of the Ohm’s law and it vanishes exactly, in 
the plasma frame, when the conductivity goes formally to inhnity. In the same limit the 
contribution of the chemical potential to the anomalous hypermagnetic diffusivity equation 
gets always erased. At hnite conductivity the anomalous contribution can be often rephrased 
in terms of the magnetic gyrotropy [ 2 ] which dehnes hypermagnetic knot solutions [TO] . 

B.l Hypermagnetic knots 

The conhgurations minimizing the hypermagnetic energy density with the constraint that 
the helicity be conserved coincide, in the perfectly conducting limit, with the ones obtainable 
in ideal magnetohydrodynamics where the anomalous currents are neglected [Dial 13111]. 

In the perfectly conducting limit Eq. fl5.6p leads to dtB = V x (u x 5) + 0{jij^/a) which 
is qualitatively similar to the result of conventional magnetohydrdynamics. Dehning the 
vector potential in the Coulomb gauge, the magnetic diffusivity equation becomes, up to 
small corrections, dtA = u x (V x A). The analysis of Ref. [22] can then be exploited. The 
magnetic energy density shall then be minimized in a hnite volume under the assumption 
of constant magnetic helicity by introducing the Lagrange multiplier ps- By taking the 
functional variation of ^ = /y (i%{|V x Ap — psA ■ (V x A)}, with respect to A and by 
requiring 5Q = 0, the conhgurations extremizing Q are such that VxB = psB. In performing 
the functional variation we assumed that V is the hducial volume of a closed system. 

The conhgurations V x B = psB have been used to describe hypermagnetic knots (see 
[lOj . third and fourth papers); in this case q has dimensions of an inverse length and sets the 


cos 6 * 
2 rA 
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2£|| 


+ ne‘i{eL-eR)sm.9 + n 

2^4 


- £r) + £4{£l + £r) - cos^ 9 
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scale of the hypermagnetic knot which is related to Chern-Simons waves. The conhgurations 
with constant pb represent the lowest state of magnetic energy which a closed system may 
attain also in the case where anomalous currents are present, provided the ambient plasma 
is perfectly conducting. 

B.2 Gyrotropic bases 

The knotted solutions can be expanded in an appropriate gyrotropic basis. Let us then 
consider a vector held a helds satisfying a x (V x a) = 0. The simplest realization of 
these Beltrami helds is provided by the eigenvectors of the curl operator but more general 
situations are know and have been extensively examined in the literature. Two gryrotropic 
and orthonormal bases of opposite parity are given by (a, 6, z) and by (c, d, z) 

d{z,p) = {cospz, — sinpz,0}, b{z,p) = {sinp^;, cospz,0}, (B.l) 

c{z,p) = {cospz, smpz,0}, d{z,p) = {—sinpz, cospz, 0}, (B.2) 

As anticipated the bases of Eqs. (IB.lj) and flB.2j) are orthonormal. Indeed we have d ■ b = 
d ■ z = b ■ z = 0 and {d x b) ■ z = 1 (and similarly for c, d and z). The unit vectors of Eqs. 
flB.ip and flB.2p are normalized eigenvectors of the curl operator with eigenvalues +p and 

-p. 

In ordinary MHD knotted solutions can be constructed from Beltrami helds by postulat¬ 
ing a solenoidal (static) current and by neglecting the displacement current. In anomalous 
magnetohydrodynamics these simple constructions cannot be immediately extended because 
of the magnetic and vortical currents. The knot solutions obtainable by extremizing the 
functional Q correspond to uniform magnetic helds well inside the core of the knot. This 
conclusion is evident if we use the basis of Eqs. (IB.ip and flB.2p . For instance in the limit 
pz < 1 the held conhguration B{z,p) = Bod{z,p) Bqx is practically uniform and directed 
along the x axis. The connections between Beltrami helds, force-free solutions in ordinary 
MHD equilibrium and electromagnetic waves propagation have been explored in a number of 
papers [22l [23l 1201E] starting from the classic works of Fermi and Chandrasekhar [28] . It is 
also possible to obtain hypermagnetic knot solutions with hnite helicity and hnite gyrotropy 
which do not satisfy the relation of Beltrami helds. These solutions have been studied in a 
number of interesting frameworks (see last two papers of Ref. in and also [221 ED] ) ■ 
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